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We present a measurement of the charmless semileptonic decay B° — > tt~£ + is using a data sample 
containing 657 xlO 6 BB events collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider operating near the T(4S) resonance. We determine the total branching fraction of 
the decay, B(B° -> n~£+v) = (1.49 ± 0.04(stat) ± 0.07(syst)) x 10" 4 . We also report a new precise 
measurement of the differential decay rate, and extract the Cabibbo-Kobayashi-Maskawa matrix 
element using model-independent and model-dependent approaches. From a simultaneous fit 

to the measured differential decay rate and lattice QCD results, we obtain \ V u b\ = (3.43±0.33) x 1CP 3 , 
where the error includes both experimental and theoretical uncertainties. 

PACS numbers: 12.15.Hh, 13.20.He, 12.38.Qk 



Weak transitions among quark flavors in the standard 
model (SM) are described by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [l|, in which \V u b\ is one of the 
least known elements. Precise measurements of the val- 
ues of the CKM matrix elements are necessary to probe 
the quark mixing mechanism of the SM and to search for 
possible physics beyond the SM. The magnitude of the 
CKM element V u b can be determined from exclusive b — > 
ulv semileptonic decays, of which B° —> 7r~£ + v Q yields 
the most precise value for \V u b\. The differential rate of 
this decay can be expressed in terms of \V u b\ and the form 
factor f+(q 2 ), where q 2 is the square of the momentum 
transferred from the B meson to the outgoing leptons, 
q 2 = (pe + Pu) 2 [3]. The present theoretical understand- 
ing of f+(q 2 ) is limited, which is a significant source for 
systematic uncertainty in the extraction of | V u b | from this 
decay. Predictions have been obtained in unquenched lat- 
tice QCD 0,11, in light cone sum rule (LCSR) theory [|| 
and in relativistic quark models [?]|. However, these pre- 
dictions typically assume a specific shape for f + (q 2 ) and 
provide reliable predictions only in a limited q 2 range 
(lattice QCD is valid near q 2 maximum, while LCSR 
is reliable near the minimum value of q 2 ). Recently, it 
has been shown that a determination of \V u b\ indepen- 
dent of a form factor shape calculation can be achieved 
by simultaneously fitting the measured q 2 spectrum and 
lattice QCD results computed near the zero recoil of q 2 



range 0, H, resulting in \V ub \ = (3.38 ± 0.36) x 1CT 3 
using the experimental data in Ref. [Tcj for the decay 
B° — > tt~£ + i/ where the error includes both theoretical 
and experimental uncertainties. The experimental un- 
certainty is a 6% while the theoretical contribution is 
estimated to be an 8.5% Q. In addition, Ref. 11 1 re- 
ports \V u b\ = (2.95 ± 0.31) x 10~ 3 by combining mea- 
surements of B° — > it~ l + v and B + — > n a t^v; here the 
error contains a 3% contribution from the branching frac- 
tion measurement, a 5% from the shape of the q 2 spec- 
trum measured in data, and an 8.5% from the theoreti- 
cal normalization. Here we describe a study of the decay 
B° — > 7r~ l + v and measure the branching fraction and the 
q 2 spectrum. We then compare with other recent studies 



of this decay l3-15|. The differential branching frac- 
tion is measured in 13 bins of q 2 , and \V u b\ is determined 
using both model-independent and model-dependent ap- 
proaches. 



The Belle detector [16|, ll7( is a large-solid-angle mag- 
netic spectrometer that consists of a silicon vertex detec- 
tor (SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
composed of CsI(Tl) crystals (ECL) located inside a su- 
perconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the coil 
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is instrumented with resistive plate chambers to detect 
K\ mesons and to identify muons (KLM). 

The data sample corresponds to an integrated luminos- 
ity of 605 fb _1 taken at a center-of-mass (cm.) energy 
near the T(45) resonance, containing 657xl0 6 BB pairs. 
For the first sample of 152 xlO 6 BB events, an inner de- 
tector configuration with a 2.0 cm beampipe and a 3- 
layer SVD was used, while a 1.5 cm beampipe, a 4-layer 
SVD and a small-cell inner drift chamber were used to 
record the remaining 505 xlO 6 BB pairs [18j |. Another 
68 fb _1 data sample taken at a cm. energy 60 MeV be- 
low the resonance is used to study the continuum back- 
ground, e + e~ — > qq, where q — u,d, s, c. Monte Carlo 
(MC) [3, simulated events equivalent to at least ten 
times the integrated luminosity were generated to model 
the signal. Samples equivalent to ten times and six times 
the integrated luminosity were generated to simulate the 
two largest background components, b — > c decays and 
continuum, respectively. To simulate rare b — > u decays, 
samples equivalent to twenty times the integrated lumi- 
nosity were generated. Final state radiation (FSR) from 
charged particles in the final state is modeled using the 
PHOTOS package 

The decay B° — > n~ t + v is reconstructed from pairs 
of oppositely charged leptons and pions. Electron candi- 
dates are identified using the ratio of the energy detected 
in the ECL to the track momentum, the ECL shower 
shape, position matching between the track and ECL 
cluster, the energy loss in the CDC, and the response 
of the ACC counters ■ Bremsstrahlung photons emit- 
ted close to the electron direction are reconstructed and 
used to correct the electron momentum [23|]. Muons are 
identified based on their penetration range and transverse 
scattering in the KLM detector [i^. In the momentum 
region relevant to this analysis, charged leptons are iden- 
tified with an efficiency of about 90% while the probabil- 
ity to misidentify a pion as an electron (muon) is 0.25% 
(1.4%). Pion candidates are selected with an efficiency 
of 85% and a kaon misidentification probability of 19%, 
based on the responses of the CDC, ACC and TOF sub- 
detectors. All charged particles are required to originate 
from the interaction point (IP) and to have associated 
hits in the SVD. The pion and lepton candidates are fit- 
ted to a common vertex and the confidence level of the fit 
is required to be greater than 1.0%. The electron (muon) 
is required to have a laboratory frame momentum greater 
than 0.8 GeV/c (1.1 GeV/c). 

The missing energy and momentum in the cm. frame 
are defined as E miss = 2£ b cam - J2i E i an d Pmiss = 
— J2iPi> respectively, where E^ eam is the beam energy in 
the cm. frame, and the sums include all charged and 
neutral particle candidates in the event. A threshold 
energy of 50 (100) MeV is required for photon candi- 
dates in the central (side) region of the ECL. The neu- 
trino 4-momentum is taken to be p v — (|p m iss|, Pmiss), 
since the determination of p m ; ss is more accurate than 



that of the missing energy. To select events compati- 
ble with the signal decay mode, we require |Qtotai| < 
3, where Qtotai is the net charge of the event, and 
-Emiss > GeV. We denote the combined system of 
the signal pion and lepton as Y. The kinematics of 
the decay constrain the cosine of the angle between 
the B and Y directions in the cm. frame, defined by 



cos 8by — (2E^ cam Ey 



i 2 B -M Y )/(2\p B \\p Y \), where 



m B refer to the mass and mo- 



ms and \ PB \ = v/^~ 
mentum of the B meson, and Ey, My, and py refer to 
the energy, mass, and momentum of the reconstructed 
Y. Background, on the other hand, is not similarly con- 
strained. In what follows we require |cos#by| < 1. 
Signal candidates are classified by their beam-energy- 

_ . | 2 , and 

energy difference, AE = Eb c 



„ - \Pn +P£+Pu\ 

- [En+Et + Ev). Can- 
didates outside of the signal region, defined by the re- 
quirements Mbc > 5.19 GeV/c 2 and \AE\ < 1 GeV, 
are rejected. To suppress background from the contin- 
uum, the ratio of second to zeroth Fox- Wolfram mo- 
ments [2^| is required to be less than 0.35. Background 
from J/ip — > decays with one muon misidentified 

as a pion is rejected by vetoing events with a Y mass 
between 3.07 GeV/c 2 and 3.13 GeV/c 2 . The sample of 
signal candidates is divided into 13 bins of q 2 from to 

0<q 2 <16GeV 2 /c 2 




5.2 5.22 5.24 5.26 5.28 



5.22 5.24 5.26 5.28 



(GeV/c ) 



M b0 (GeV/c") 



FIG. 1: Fit projections (a,b) in AE with M bc > 5.27 GeV/c 2 , 
and (c, d) in M hc with \AE\ < 0.125 GeV. The projec- 
tions (a,c) and (b,d) show the regions q 2 < 16 GeV 2 /c 2 
and q 2 > 16 GeV 2 /c 2 , respectively. The points with er- 
ror bars are T(4S') data, the histograms are (from top to 
bottom) B° -> n~£ + 

v signal (open), B —¥ X u lv (cross- 
hatched), B — > X c lv (hatched) and continuum background 
(black-filled). The smaller error bars are statistical only while 
the larger ones include systematic uncertainties. 
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26.4 GeV 2 /c 2 (the bin width is 2 GeV 2 /c 2 , except for 
the last bin). The value of q 2 is calculated as the square 
of the difference between the 4-momenta of the B meson 
and that of the pion. As the B direction is only kinemati- 
cally constrained to lie on a cone around the Y direction, 
we take a weighted average over four different possible 
configurations of the B direction [26| . Background is fur- 
ther suppressed by applying selection criteria as a func- 
tion of q 2 to the following quantities: the angle between 
the thrust axis of the Y system and the thrust axis of 
the rest of the event; the angle of the missing momentum 
with respect to the beam axis; the helicity angle of the 
iv system [27j |: and the missing mass squared of the event, 
M mi SS = ^miss "Pmiss- The helicity angle is the angle be- 
tween the lepton direction and the direction opposite to 
the B meson in the iv rest frame. These selections are 
optimized separately in each bin of q 2 by maximizing the 
figure-of-merit S/ y/(S + B), where S (B) is the expected 
number of signal (background) events. 

The fraction of events that have multiple candidates 
is 66%. To remove multiple signal candidates in a single 
event, the candidate with the smallest iv helicity angle is 
selected. After imposing all selections described above, 
the reconstruction efficiency for signal ranges from 7.7% 
to 15.0% over the entire q 2 range. The fraction of the 
self-cross-feed component, in which one or more of the 
signal tracks are not correctly reconstructed, is 3.5%. 

The signal yield is determined by performing a two- 
dimensional, binned maximum likelihood fit to the 
(Mbc AE) plane in 13 bins of q 2 [HI. Background con- 
tributions from b — >• uiv, b —> civ and non-BB con- 
tinuum are considered in the fit. Probability density 
functions (PDFs) corresponding to these fit components 
are obtained from MC simulations. To reduce the num- 
ber of free parameters, the q 2 bins of the background 
components are grouped into coarser bins: four bins for 
b — > uiv, and three bins for b — > civ. The choice of the 
binning was chosen from the total statistical error, num- 
ber of parameters to fit, and the complexity of the fits. 
The q 2 distribution of the continuum MC [29( simulation 
is reweighted to match the corresponding distribution in 
off-resonance data. For this procedure, a continuum MC 
sample about 60 times the integrated luminosity of the 
off-resonance data is used. The continuum normaliza- 
tion is fixed to the scaled number of off-resonance events, 
52928 events. Including signal yields in each q 2 bin, there 
are 20 free parameters in the fit. 

We obtain 21486 ± 548 signal events, 52543 ± 1148 
b — > uiv events, and 161829 ± 976 b — > civ background 
events. These yields agree well with the expectations 
from MC simulation studies. The x 2 /n.d.f. of the fit is 
2962/3308. The projections of the fit result in AE and 
Mbc are shown in Fig. [T] for the regions q 2 < 16 GeV 2 /c 2 
and q 2 > 16 GeV 2 /c 2 . Bin-to-bin migrations due to 
q 2 resolution are corrected by applying the inverse detec- 
tor response matrix [30l ] to the measured partial yields. 



The partial branching fractions AB are calculated us- 
ing the signal efficiencies obtained from MC simulation. 
The total branching fraction B is the sum of partial 
branching fractions taking into account correlations when 
calculating the errors. We find B(B° — > ir~i + v) = 
(1.49±0.04(stat)±0.07(syst)) x 10~ 4 , where the first er- 
ror is statistical and the second error is systematic. This 
result is significantly more precise than our previous mea- 
surement jlH with B -> D^t+v tags on a 253 fb" 1 data 
sample. 

To estimate the systematic uncertainties on AB, we 
include the following contributions: the uncertainties in 
lepton and pion identification, the charged particle re- 
construction, the photon detection efficiency, and the re- 
quirement on the x 2 probability of the vertex fit, which 
is estimated by comparing results with and without this 
requirement. The results are summarized as detector ef- 
fects in Table [I] They depend weakly on q 2 and amount 
to 3.4% for the entire q 2 range. We vary the branching 
fractions of the decays contributing to the b —y uiv and 
b — » civ backgrounds within ±1 standard deviation of 
their world-average values [31| and assign an uncertainty 
of 0.6% to the total yield. We further consider form fac- 
tor uncertainties in the decays B° — > ir~i + v 
-i+v BH, B a -> D-i+v and B° -> D* 



14, 



B° -> 

p~l^v [g, ±S" -> D~l^v and tS" -» Lr~i + v [33| . 
and uncertainties in the shape function parameters of 
the inclusive b —5- uiv model [341. These uncertainties 



x10 




15 20 25 
Unfolded q 2 (GeV 2 /c 2 ) 



FIG. 2: Distribution of the partial branching fraction as 
a function of q 2 after unfolding (closed circles). The er- 
ror bars show the statistical and the total uncertainty on 
the data. The curve is the result of a fit of the BK form 
factor parameterization [35j to our data. The four his- 
tograms (dashed:ISGW2; plaimHPQCD; dotted:FNAL; dot- 
dashed:LCSR) show various form factor predictions. 
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TABLE I: Values of AB(q 2 ) and relative uncertainties (%). The uncertainties in MC input parameters are given separately for 
branching fractions (BF) and form factors (FF). 



g 2 (GeV 2 /c 2 ) 


- 6 


6 - 12 


12 - 18 


18 - 26.4 


- 16 


16 - 26.4 


Total 


AB (x 10 7 ) 


391.19 


434.25 


389.47 


279.18 


1096.34 


397.75 


1494.09 


Detector effects 


3.4 


3.5 


3.5 


3.5 


3.4 


3.5 


3.4 


Physics parameters (BF) 


0.8 


0.7 


0.6 


0.7 


0.6 


0.6 


0.6 


Physics parameters (FF) 


1.9 


1.7 


1.9 


1.8 


1.3 


1.8 


1.1 


Continuum correction 


4.4 


2.3 


3.4 


2.3 


2.1 


2.6 


1.8 


Other sources 


2.1 


2.5 


2.4 


2.4 


2.1 


2.3 


2.0 


Total statistical error 


5.3 


3.9 


4.8 


6.1 


3.0 


5.3 


2.6 


Total error 


8.2 


6.5 


7.5 


8.1 


5.7 


7.5 


5.2 



correspond to a 1.1% error on B(B° — > n~l + v). The 
uncertainty in the correction of the continuum MC is 
estimated by varying its weights by their statistical un- 
certainties. The other sources of systematic uncertainty 
in Table U include the uncertainty in the T(45) -> B a B a 
branching fraction [3l|, limited MC statistics, the effect 
of final state radiation, which is estimated by investigat- 
ing MC samples with and without bremsstrahhmg cor- 
rections calculated using the PHOTOS package, and the 
uncertainty in the number of BB pairs in the data sam- 
ple. For values of AB in individual q 2 bins, a breakdown 
of the systematic uncertainties is presented in Table [V] 
and the statistical and systematic correlations is given in 
Table HI] and Table [IVJ 

We fit the AB distribution using the two-parameter 
BK parameterization [35j of f+(q 2 ), taking into account 
statistical and systematic correlations. The result is 
shown in Fig. [2l Although this parameterization has 
been criticized j36| , we present the fit result in order to di- 
rectly compare with other existing results [l(| • We obtain 
\V ub \f + {0) = (9.24±0.18(stat) ± 0.21(syst)) x 10" 4 and 
a = 0.60 ± 0.03(stat) ± 0.02(syst), where a is a positive 
constant that scales with tub [351 ] . The x 2 probability 
of the fit is 62%. We also calculate the % 2 probabilities 
of different theoretical form factor predictions with our 
binned data. We obtain probabilities of 42% and 43% 
for the HPQCD [4] and the FNAL 0] lattice QCD cal- 
culations, respectively, and 49% for the LCSR theory @. 
The ISGW2 quark model . 7J, for which the probability is 
2.3 xl0~ 6 , is incompatible with the experimental data. 

As described in Ref. [9(, the CKM matrix element \V u b\ 
can be extracted from a simultaneous fit to experimental 
and lattice QCD results (from the FNAL/MILC Collabo- 
ration 0), taking into account statistical and systematic 
correlations. To this end, the o 2 variable is transformed 
to a dimensionless variable z [a . In addition, the two 
functions, P+ and 4>+ are taken from Ref. 37], where P. 



+ 



b + are taken from Ref. 
is a function that accounts for the pole at q 2 = m 2 B , and 
<fi + is an analytic function that controls the values of the 
on series coefficients. In terms of the new variable z, the 



product of the form factor f+(q 2 ) and the functions P+ 
and 4> + has the simple form, Y^Lq QiZ ' ''■ We n * the lattice 
QCD results and experimental data with a third-order 
polynomial where the free parameters of the fit are the 
coefficients at and the relative normalization between lat- 
tice QCD results and experimental results, which is \V u b\- 
The resulting experimental data (which are scaled by the 
fitted \ V u b\ value) and the lattice QCD results are shown 
in Fig. [3] We obtain |V u6 | = (3.43 ± 0.33) x 10~ 3 , a = 
0.022 ±0.002, a x = -0.032 ± 0.004, a 2 = -0.080 ±0.020 
and a 3 = 0.081 ± 0.066, where the x 2 /n.d.f. of the fit is 
approximately 12/20. Statistically, we find no significant 
difference in the fitted value of \V u b\ using second- and 
fourth-order polynomial fits. Note that the error in \ V u b\ 
includes both experimental and theoretical uncertainties. 
We find that the error includes a 3% contribution from 
the branching fraction measurement, a 4% from q 2 shape 



0.04 r 



> 0.035 1 



-0 — MILC lattice points 

_0 Belle data scaled by V J from fit 

Simultaneous fit 




-0.3 -0.2 -0.1 



FIG. 3: \V u b\ extraction from a simultaneous fit of exper- 
imental (closed circles) and FNAL/MILC lattice QCD re- 
sults (open circles) [jj. The error for each experimental data 
point is the total experimental uncertainty. The smaller error 
bars of the lattice QCD results are statistical only while the 
larger ones also include systematic uncertainties. 
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TABLE II: Values extracted for \V u b\ using different form fac- 
tor predictions. The first error on \V u b\ is the experimental 
error including statistical, systematic uncertainties and the 
uncertainty in the B° lifetime [U, the last asymmetric er- 
rors arise from the uncertainty in A£. 





q 2 (GeV 2 /c 2 ) 


AC (ps- 1 ) 


\V ub \ (10- 3 ) 


HPQCD [4] 


> 16 


2.07 ±0.57 


3.55 ± 0.13±°;g 


FNAL [5] 


> 16 


1.83 ±0.50 


3.78 ± 0.14±g;|| 


LCSR [6J 


< 16 


5.44 ± 1.43 


3.64±0.1lloio 



measured in data, and an 8% uncertainty from theoreti- 
cal normalization. The experimental and the total errors 
are compatible with the previous results in Ref. 0, [ll[ • 

Alternatively, \V u b\ can be determined from the mea- 
sured partial branching fraction using the relation \V u b \ = 
^AB/(t b oA(), where tb« is the B° lifetime [3l( and 
AC is the normalized partial decay width derived in dif- 
ferent theoretical approaches [4j46|. These calculations 
typically assume a specific parameterization of the form 
factor shape. Values of \V u b\ for different form factor pre- 
dictions are given in Table HU 

In summary, using 657xl0 6 BB events of Belle 
T(4S') data we measure the partial branching fractions 
of the decay B° — > tt~1 + v in 13 bins of q 2 . The to- 
tal branching fraction is found to be (1.49 ± 0.04(stat) ± 
0.07(syst)) x 10~ 4 . A combined fit of experimental and 
FNAL/MILC lattice QCD results @, yields a new pre- 
cise determination of \V u b\ from this decay, \V u b\ = 
(3.43 ± 0.33) x 10~ 3 . Determinations using only a frac- 
tion of the phase space lead to less precise but sta- 
tistically compatible numbers for \V u b\- using a LCSR 
calculation for the region q 2 < 16 GeV 2 /c 2 Q yields 
(3.64±0.06(stat)±0.09(syst)±gjg(FF)) x 1CT 3 . Assum- 
ing the HPQCD @ and the FNAL 0] lattice QCD cal- 
culations, sensitive to the region q 2 > 16 GeV 2 /c 2 , we 
obtain (3.55 ± 0.09(stat) ± 0.09(syst)t^ 2 (FF)) x 10~ 3 
and (3.78 ± O.lO(stat) ± 0.10(syst)t£g(FF)) x lO" 3 , re- 
spectively. 
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TABLE III: Correlation coefficients of the AB(B° — > ty £ + is) statistical errors. 
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TABLE IV: Correlation coefficients of the AB(B° — ¥ it £ + v) systematic errors. 
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TABLE V: The raw, unfolded yields, signal efficiencies and the partial branching fractions and their relative errors (%) from various sources in 13 bins of q 2 . 
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